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Abstract. We study the spectra and pseudospectra of semi-infinite and bi-infinite tridiagonal random 
matrices and their finite principal submatrices, in the case where each of the three diagonals varies 
over a separate compact set, say U,V,W C C. Such matrices are sometimes termed stochastic Toeplitz 
matrices A+ in the semi-infinite case and stochastic Laurent matrices A in the bi-infinite case. Their 
spectra, E = spec A and E+ = specA+, are independent of A and A+ as long as A and A+ are 
pseudoergodic (in the sense of E.B. Davies, Commun. Math. Phys. 216 (2001), 687-704), which holds 
almost surely in the random case. This was shown in Davies (2001) for A; that the same holds for A^ 
is one main result of this paper. Although the computation of E and E 4 . in terms of U, V and W is 
intrinsically difficult, we give upper and lower spectral bounds, and we explicitly compute a set G that 
fills the gap between E and E+ in the sense that E U G = E+. We also show that the invertibility of 
one (and hence all) operators A+ implies the invertibility - and uniform boundedness of the inverses 
- of all finite tridiagonal square matrices with diagonals varying over U, V and W. This implies that 
the so-called finite section method for the approximate solution of a system Aj^x = b is applicable 
as soon as A+ is invertible, and that the finite section method for estimating the spectrum of A+ 
does not suffer from spectral pollution. Both results illustrate that tridiagonal stochastic Toeplitz 
operators share important properties of (classical) Toeplitz operators. Indeed, one of our main tools is 
a new stochastic version of the Coburn lemma for classical Toeplitz operators, saying that a stochastic 
tridiagonal Toeplitz operator, if Fredholm, is always injective or surjective. In the final part of the 
paper we bound and compare the norms, and the norms of inverses, of bi-infinite, semi-infinite and 
finite tridiagonal matrices over U, V and W. This, in particular, allows the study of the resolvent 
norms, and hence the pseudospectra, of these operators and matrices. 
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Keywords: finite section method, random operator, Coburn lemma, Jacobi operator, spectral 
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1 Introduction and Main Results 


In this paper we study so-called Jacobi operators over three sets U, V and W, meaning bi- and 
semi-infinite matrices of the form 
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with entries ut € U, Vi £ V and Wi € W for all i under consideration. The sets U, V and W are 
nonempty and compact subsets of the complex plane C, and the box marks the matrix entry of A at 
(0, 0). We will be especially interested in the case where the matrix entries are random (say i.i.d.) 
samples from U, V and W. Trefethen et al. [57] call the operator A a stochastic Laurent matrix 
in this case and A^ a stochastic Toepiitz matrix. We will adopt this terminology which seems 
appropriate given that one of our aims is to highlight parallels between the analysis of standard 
and stochastic Laurent and Toepiitz matrices. 

It is known that the spectrum of A depends only on the sets U, V, and W, as long as A is 
pseudoergodic in the sense of Davies [20], which holds almost surely if A is stochastic (see the 
discussion below). Via a version, which applies to stochastic Toepiitz matrices, of the famous 
Coburn lemma [18] for (standard) Toepiitz matrices, a main result of this paper is to show that, 
with the same assumption of pseudoergodicity implied by stochasticity, also the spectrum of A+ 
depends only on [/, V, and W. Moreoever, we tease out very explicitly what the difference 
is between the spectrum of a stochastic Laurent matrix and the spectrum of the corresponding 
stochastic Toepiitz matrix. (The difference will be that certain ‘holes’ in the spectrum of the 
stochastic Laurent case may be ‘filled in’ in the stochastic Toepiitz case, rather similar to the 
standard Laurent and Toepiitz cases.) 

A second main result is to show that infinite linear systems, in which the matrix, taking one 
of the forms (1), is a stochastic Laurent or Toepiitz matrix, can be solved effectively by the 
standard finite section method, provided only that the respective infinite matrices are invertible. 
In particular, our results show that, if the stochastic Toepiitz matrix is invertible, then every finite 
n X n matrix formed by taking the first n rows and columns of A+ is invertible, and moreover 
the inverses are uniformly bounded. Again, this result, which can be interpreted as showing that 
the finite section method for stochastic Toepiitz matrices does not suffer from spectral pollution 
(cf. [42]), is reminiscent of the standard Toepiitz case. 

Related work. The study of random Jacobi operators and their spectra has one of its main 
roots in the famous Anderson model [1, 2] from the late 1950’s. In the 1990’s the study of a non- 
selfadjoint (NSA) Anderson model, the Hatano-Nelson model [31, 44], led to a series of papers on 
NSA random operators and their spectra, see e.g. [25, 20, 19, 43]. Other examples of NSA models 
are discussed in [16, 17, 57, 37]: one example that has attracted significant recent attention (and 
which, arguably, has a particularly intriguing spectrum) is the randomly hopping particle model 
due to Feinberg and Zee [21, 22, 11, 13, 12, 27, 28, 29, 15]. A comprehensive discussion of this 
history, its main contributors, and many more references can be found in Sections 36 and 37 of 
[58]. 

A theme of many of these studies [57, 58, 11, 13, 12, 28, 15], a theme that is central to this 
paper, is the relationships between the spectra, norms of inverses, and pseudospectra of random 
operators, and the corresponding properties of the random matrices that are their finite sections. 
Strongly related to this (see the discussion in the ‘Main Results’ paragraphs below) is work on the 
relation between norms of inverses and pseudospectra of finite and infinite classical Toepiitz and 
Laurent matrices [51, 3, 5]. In between the classical and stochastic Toepiitz cases, the same issues 
have also been studied for randomly perturbed Toepiitz and Laurent operators [8, 9, 10, 7]. This 
paper, while focussed on the specific features of the random case, draws strongly on results on the 
finite section method in much more general contexts: see [35] and the ‘Finite sections’ discussion 
below. With no assumption of randomness, the finite section method for a particular class of NSA 
perturbations of (selfadjoint) Jacobi matrices is analysed recently in [42]. 

We will build particularly on two recent studies of random Jacobi matrices A and A+ and their 
finite sections. In [27] it is shown that the closure of the numerical range of these operators is the 
convex hull of the spectrum, this holding whether or not A and A_|_ are normal operators. Further, 
an explicit expression for this numerical range is given: see (33) below. In [39] progress is made 
in bounding the spectrum and understanding the hnite section method applied to solving infinite 
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linear systems where the matrix is a tridiagonal stochastic Laurent or Toeplitz matrix. This last 
paper is the main starting point for this present work and we recall key notations and concepts 
that we will build on from [39] in the following paragraphs. 

Matrix notations. We understand A and A+ as linear operators, again denoted by A and A+, 
acting boundedly, by matrix-vector multiplication, on the standard spaces and £^(N) of bi- 

and singly-infinite complex sequences with p S [1, oo]. The sets of all operators A and A^ from (1) 
with entries Ui G U, Vi G V and Wi G W for all indices i that occur will be denoted by M{U, V, W) 
and M+{U, V^W), respectively. For n £ N the set of n x n tridiagonal matrices with subdiagonal 
entries in U, main diagonal entries in V and superdiagonal entries in W (and all other entries 
zero) will be denoted Mn{U, V, IT), and we set MuniU, V, W) := V, W). For X = £^(1) 

with I = N or Z, we call a bounded linear operator A : X ^ X a band-dominated operator and 
write A G BDO(X) if A is the limit, in the operator norm on X, of a sequence of band operators 
(i.e. bounded operators on X that are induced by infinite matrices with finitely many nonzero 
diagonals). For A = (aij)ijgi G BDO(X), A^ = {aji)ij^i G BDO(X) denotes the transpose 
of A. The boundedness of the sets U, V and W implies that every operator in M{U,V,W) or 
M^{U, V, W) is bounded and hence band-dominated (of course even banded). We use || • || as the 
notation for the norm of an element of X = f^(I), whether I = N, Z, or —N = {..., —2, —1}, or I is 
a finite set, and use the same notation for the induced operator norm of a bounded linear operator 
on X (which is the induced norm of a finite square matrix if I is finite). If we have a particular 
p G [1, oo] in mind, or want to emphasise the dependence on p, we will write jj • jjp instead of jj ■ jj. 

Random alias pseudoergodic operators. Our particular interest is random operators in 
M{U, V, W) and M+{U, V, W). We model randomness by the following deterministic concept: we 
call A G M{U,V,W) pseudoergodic and write A G 'I'E([/, V, W) if every finite Jacobi matrix over 
U, V and W can be found, up to arbitrary precision, as a submatrix of A. Precisely, for every B G 
V, W) and every e > 0 there is a finite square submatrix C oi A such that ||R —Cj] < er. Here 
C is composed of rows i = /c-l-l,..., k-\-n and columns j = ..., k-\-n of A, where k is some integer 

and n is the size of B. By literally the same definition we define semi-infinite pseudoergodic matrices 
and denote the set of these matrices by 'I'E+(t/, V, IF). Pseudoergodicity was introduced 
by Davies [20] to study spectral properties of random operators while eliminating probabilistic 
arguments. Indeed, if all matrix entries in (1) are chosen independently (or at least not fully 
correlated) using probability measures whose supports are U, V and IF, then, with probability 
one, A and A^ in (1) are pseudoergodic. 

Fredholm operators, spectra, and pseudospectra. Recall that a bounded linear operator 
B : X ^ Y between Banach spaces is a Fredholm operator if the dimension, a{B), of its null- 
space is finite and the codimension, P{B), of its image in Y is finite. In this case, the image of 
B is closed in Y and the integer indS := a{B) — fl(B) is called the index of B. Equivalently, 
i? is a Fredholm operator if it has a so-called regularizer C : F —X modulo compact operators, 
meaning that BC — ly and CB — Ix are both compact. For a bounded linear operator B on 
£^(1) with I G {Z,N, —N}, we write spec^H and specf^g i? for the sets of all A £ C for which 
B — XI is, respectively, not invertible or not a Fredholm operator on £p(I). Because A and 
in (1) are band matrices, their spectrum and essential spectrum do not depend on the underlying 
^^’-space [34, 36, 52] , so that we will just write spec A and speCggg A for operators A G M{U, V, IF) 
- and similarly for A+ G M+(U, V, IF). Even more, if A+ G M+(U^ V, IF) or H £ M{U, F, IF) is 
Fredholm as an operator on £p(N) or £p(Z), respectively, both its Fredholm property and index are 
independent oi p [52, 36]. 

Following notably [58], we will be interested in not just the spectrum and essential spectrum, 
but also the e-pseudospectrum, the union of the spectrum with those A £ C where the resolvent 
is well-defined with norm > e~^. Precisely, for an n x n complex matrix B, or a bounded linear 
operator B on £P{Tj with I £ {Z,N}, we define, for I < p < oo and e > 0, 

spec^H := {A £ C : \\{B - XI)-^\\p > e-^} , (2) 
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with the convention that ||A“^||p := oo if A is not invertible (so that spec^i? C specfi?). While 
specks is independent of p (and so abbreviated as spec B), the set specfi? depends on p in general. 
It is a standard result (see [58] for this and the other standard results we quote) that 

spec i? + eD c speeds, (3) 

with D := {z G C : jzj < 1} the open unit disk. If p = 2 and B is & normal matrix or operator 
then equality holds in (3). Clearly, for 0 < ei < £ 2 , specB C spec^^B C spec^^i?, and speci? = 
n£>oSpecfi3 (for 1 < p < 00 ). Where S denotes the closure of S' C C, a deeper result, see the 
discussion in [54] (summarised in [12]), is that 

speefS = Speeds := {A G C : 1](B - A/)-il|p > £-^} . (4) 

Interest in pseudospectra has many motivations [58]. One is that speefi? is the union of spec(i?+T) 
over all perturbations T with JjTjjp < £. Another is that, unlike speci? in general, the pseudospec¬ 
trum depends continuously on B with respect to the standard Hausdorff metric (see (63) below). 

Limit operators. A main tool of our paper, and of [39], is the notion of limit operators. 
For A = G BDO(A') with X = P'{'L) and in Z with —>■ 00 we say that 

B = {bij)ij^z is a limit operator of A if, for all i,j G Z, 

ai+h„,j+h^ ^ bij as n-)• 00 . (5) 

The boundedness of the diagonals of A ensures (by Bolzano-Weierstrass) the existence of such 
sequences {hn) and the corresponding limit operators B. From A G BDO(A) it follows that 
B G BDO(X). The closedness of U, V and W implies that B G M{U,V,W) if A G M{U,V,W). 
We write (t°p(A) for the set of all limit operators of A. Similarly, B = {bij)ij^’z is a limit operator 
of A_|_ = (ay )ijgN £ BDO(£P(N)) if (5) holds for a sequence (/i„) in N with hn -too. Note that 
limit operators are always given by a bi-infinite matrix, no matter if the matrix A or A+ to start 
with is bi- or semi-infinite. The following lemma summarises the main results on limit operators: 

Lemma 1.1 Let A G BDO(^^(I)) with I G {Z, N} and let B be a limit operator of A. Then: 

a) [ 47 ] It holds that jjBjj < JJA]]. 

b) [ 47 , 49 ] If A is Fredholm then B is invertible, and B~^ is a limit operator of any regularizer 
C of A. (Note that B,B-^ G BDO(^p(Z)) and C G BDO(^p(I)) hold if A & BDO(F(I)).; 

c) [ 47 , 35, 14 , 40 ] A is Fredholm iff all its limit operators are invertible. 

d) [20, 35] If A G M{U,V,W) or A G Mj^{U,V,W) then A is pseudoergodic iff a°'^{A) = 
M{U,V,W). 

e) [ 46 ] If A G BDO(£^(N)) is Fredholm then B+ := {Bij)ij^fi is Fredholm and ind (A) = 
ind (B+). 

So we immediately get that A G 'I'E([/, V, W) is Fredholm iff all i? G M(C/, V, W) are invertible, in 
which case of course A G M{U, V, W) is invertible. 

Finite sections. A further topic of [39] and our paper is the so-called finite section method 
(FSM). This method aims to approximately solve an equation Ax = b, i.e. 

x{j) = bii), iGZ, (6) 

iez 

by truncating it to 

^ ^ C’ij — ^(^5 i rnj (f) 

ln.<0<rn 
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where the cut-off points Zi,Z 2 ,--- —oo and ri,r 2 ,... —>■ -|-oo are certain, sometimes well-chosen, 
integers. The aim is that, assuming invertibility of A (i.e. unique solvability of (6) for all right- 
hand sides b), also (7) shall be uniquely solvable for all sufficiently large n and the solutions x„ 
shall approximate the solution x of (6) (in the sense that, for every right hand side 6, it holds as 
n —>■ oo that \\xn — x|| 0, if 1 < p < oo, that ||a;„|| = 0(1) and x„(j) x{j) for every j G Z, if 

p = oo). If that is the case then the FSM is said to be applicable to A (or we say that the FSM 
applies to A). 

If i and j in (6) only run over the positive integers, N, then this system corresponds to a 
semi-infinite equation A+x+ = 6+. The FSM is then to freeze at 1 and only let r„ go to -boo. 
Otherwise the terminology is identical. 

Applicability of the FSM is equivalent [45, 53, 49] to invertibility of A plus stability of the 
sequence of finite matrices 

An ■= n=l,2,... . (8) 

The latter means that, for all sufficiently large n, the matrices (the so-called finite sections of A) 
are invertible and their inverses are uniformly bounded, in short: limsup„_,,oo ||A“^|| < oo. This, 
moreover, is known [39, 55] to be equivalent to the invertibility of A and of certain semi-infinite 
matrices that are associated to A and to the cut-off sequences (/„) and (r„). Those associated 
semi-infinite matrices are partial limits (in the strong topology) of the upper left and the lower 
right corner of the finite matrix as n ^ oo. Precisely, the associated matrices are the entrywise 
limits 


5+ and C_ as n —> oo (9) 

of semi-infinite submatrices of A, where {ln)'^=i and (r(j))]T 2 are subsequences of and 

(r„)))Ti, respectively, such that the limits (9) exist. So B+ and C- are one-sided truncations of 
limit operators of A; they tell us what we find in the limit when jumping along the main diagonal 
of A via the sequences ^i, h ,... and ri,r 2 , ... - or subsequences thereof. Hence, by the choice of the 
cut-off sequences (/„) and (r„), one can control the selection of associated matrices B+ and C- 
and consequently control the applicability of the FSM. Let us summarize all that: 

Lemma 1.2 For A = € BDO(£p(Z)) and two cut-off sequences {ln)ffLi {,T'n)ff=i in 

Z with In —oo and Xn -boo, the following are equivalent: 

i) the FSM (7) is applicable to A, 
a) the sequence (An)ffLi, niith An from (8), is stable, 

Hi) A and all limits B+ and C- from (9) are invertible. 

iv) A and all limits B+ and C- from (9) are invertible, and the inverses Bf^ and CZ^ are 
uniformly bounded. 

Proof. That applicability i) is equivalent to invertibility of A plus stability ii) is a classical result 
(called “Polski’s theorem” in [26]) for the case of strong convergence A„ —>■ A, and it is in [53] 
for the more general case considered here. (Note that the convergence An —)• A is generally not 
strong if p = oo.) That ii) is equivalent to iv) was shown in [48] for p = 2, in [49] for p G [I, oo) 
and in [35] for p G [l,oo]. So in particular, ii) implies invertibility of A and hence also i). The 
equivalence of Hi) and iv) is shown in [50] for p G (I,oo) and in [38] for p G [l,oo]. The case of 
arbitrary monotonic cut-off sequences (Z„) and (r„) can be found in [55, 39]. ■ 

If the FSM is applicable to A then, for every right hand side b, Xn(j) x(j) as n —)• oo for 
every j G Z, where Xn is the solution to (7) and x that of (6). But this implies that j|a:|| < 
liminf„_^oo ll^Cnjl, for every b, and hence 

liminf and similarly liminfjjA„jj > jjAjj, (10) 
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this latter holding whether or not the FSM is applicable to A. Complementing this bound, we 
remark that it can be deduced from [56, Proposition 3.1] (or see [30, Section 6]) that 

hmsnp p-i|| = snp{p-i||, \\B^% lie'll}, (11) 

n—^oo 

where this supremum is taken over all limits B+ and C- in (9), and is attained as a maximum 
if the FSM is applicable to A. In our arguments below we will not need (11), however, only the 
much simpler lower bound (10). 

Versions of Lemma 1.2, (10), and (11) hold for semi-infinite matrices A^ = (aij)ijgN S 
BDO(£P(N)), with the modification that In = 1, which implies that every limit B+ in (9) is 
nothing but the matrix A+ again, so that in Lemma 1.2 hi), iv) and (11) it is the invertibility of 
only yl+ and C_ that is at issue. 

Remark 1.3 ~ Reflections. Often we find it convenient to rearrange/reflect the matrices C- = 
{cij)i j=_ao from (9) as 5+ = (c_j^_i)^_o- This rearrangement C- i—>■ 5+ corresponds to a matrix 
reflection against the bi-infinite antidiagonal; it can be written as B+ = RC^R, where R denotes 
the bi-infinite flip {xi)i^z As an operator on one gets ||R-i-||p = ||i?Cji?||p = 

IjCJlIp = IlC-llq with^ p~^ + q~^ = 1. When we speak below about the FSM for A and its 
“associated semi-infinite submatrices R+” we will mean all R+ from (9) plus the reflections R+ = 
RCZR of all C_ from (9). □ 

A simple choice of cut-off sequences is to take In = —n and Vn = n for n = 1,2,.... This is 
called the full FSM for A. For a semi-infinite matrix A+ the full FSM is to take In = f and rn = n 
for n = 1, 2,.... In either case, the full FSM leads to more associated matrices R+ (and hence 
to a smaller chance for applicability of the FSM) than “thinning out” those cut-off sequences 
in a way that suits the matrix A (or A+) at hand. For example, if A G 'FF{U,V,W) then all 
R+ G M_|_(?7, y, TV) are associated to A in case of the full FSM. So, in addition to A itself, all 
R+ G M+{U, V, W) have to be invertible to make sure the full FSM applies to A. That is why, 
in [39], the cut-offs In and rn have been placed very sparsely and in a special way that leads to 
all associated R+ being Toeplitz. A simple consequence of Lemma 1.2 is the following lemma 
which also (trivially, by the equivalence of i) and ii) in Lemma 1.2) holds for semi-infinite matrices 
A+ = (aij)ijgN (with = 1 in that case). 

Lemma 1.4 If the full FSM applies to A = G BDO{£^{'Z)) then the FSM with any 

monotonic cut-off sequences {ln)fliLi and (r„)(]V]^ applies to A. 

Proof. Suppose that the full FSM method is applicable to A. Fix two arbitrary monotonic cut-off 
sequences {ln)^=i —oo and {rn)^=i —t +oo and look at two associated matrices R+ and C- 
from (9) with respect to subsequences (I'n) and (r(j) of (1„) and (r„). Since (/„) and (r„) are (at 
least for sufficiently large n) subsequences of the sequences (—1, —2, —3,...) and (1,2,3,...) that 
are used for the full FSM, the same is true for {ff) and (r(j). So 5+ and C- are also associated 
to A in case of the full FSM and so, by the equivalence of i) and iii) in Lemma 1.2, are invertible 
together with A. Again by the equivalence of i) and iii) in Lemma 1.2, since A and all matrices 
R+ and C- associated to these cut-off sequences are invertible, the FSM with cut-off sequences 
(/„) and (r„) applies to A. ■ 

Lemma 1.4 is why we place particular focus on the full FSM: it is the most demanding version 
of the FSM - if this version applies then all ent-off sequences will be fine. 

Main Results. Having set the notations, let us now sketch our main results. For operators 
A G 4>E(?7, y, W) and R+ G 4/E_|_({7, V, VF), we are interested in the four sets 

speCggg A., spec A, speCgssi?+ and speci?+. 

^Here and in what follows we put oo“^ := 0, so that p = 1 =1- g = oo and p = oo => q = 1. 
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From [39] we know that the first three sets coincide, 

speCe,,A = spec^ = speCe,,B+, (12) 

and are independent of A and i3+, as long as these are pseudoergodic. We will show that also the 
fourth set, speci?+, is independent of and we indicate what the difference between the two 
sets is. The key to describe the difference between speCg^g and speci3-|- is a new result that 
has a famous cousin in the theory of Toeplitz operators: Coburn’s Lemma [18] says that, for every 
bounded and nonzero Toeplitz operator T+, one has a(7+) = 0 or /3(T'+) = 0, so that T+ is known 
to be invertible as soon as it is Fredholm and has index zero. We prove that the same statement 
holds with the Toeplitz operator r+ replaced by any € M+{U, V, W) provided that 0 is not in 
(12). So speCgssil+ and specil+ differ by the set of all A € C for which B+ — A/+ is Fredholm 
with a nonzero index. We give new upper and lower bounds on the sets speCg^g B+ and specil+, 
and we find easily computable sets G that close the gap between the two, i.e., sets G for which it 
holds that spec B+ = G U speCg^g B+. 

On the other hand, knowledge about invertibility of semi-infinite matrices i?+ € M+{U, V, W) 
is all we need to study applicability of the FSM, so that our new Coburn-type result has immediate 
consequences for the applicability of the FSM (even the full FSM) to pseudoergodic operators. In 
[39] the question of the applicability of the full FSM to an operator G '1>E+(17, V, W) could not 
be settled (nor could, in [39], the applicability of the full FSM to A G 'I'E(17, V, W)-. for brevity we 
just focus on the semi-infinite case in this paragraph). Instead, the cut-off sequences for the FSM 
were chosen (“adapted to A+”) in a way that made all associated semi-infinite matrices G_ in (9) 
Toeplitz. Classical Coburn then implied invertibility of all these G_, as a consequence of their 
being Fredholm of index zero, which holds as long as A-f- is invertible. Thanks to the new Coburn 
result, this “adaptation” twist is no longer needed. The full FSM can be seen to apply by exactly 
the same argument, but with the associated operators G_ no longer required to be Toeplitz. 

Perhaps one of the main messages of our paper is that operators in 4'E+(17, E, W) (termed 
“stochastic Toeplitz operators” in [57]) behave a lot like usual Toeplitz operators when it comes to 

• the gap between essential spectrum and spectrum (both enjoy a lemma of Coburn type), and 

• having an applicable FSM (in both cases, the FSM applies iff the operator is invertible). 

Similar coincidences can be shown for operators in 'I'E(C/, E, W) (the “stochastic Laurent operators” 
in the terminology of [57]) and usual Laurent operators. 

In Section 3 of our paper we show that the full FSM applies to A^ G 4'E_|_(C/, E, IE), and 
automatically to all other operators in M+{U, V, IE) and in M{U, E, IE), as soon as A^ is invertible. 
Even more, we show that all matrices in Mfin{U,V,W) are invertible if A+ is invertible, so that 
the truncated systems (7) are uniquely solvable for all n > I (as opposed to n > no with an 
no that nobody knows). If A^ G 'I'E+(C/, E, IE) is not invertible but Fredholm with index k = 
k(C/, E, IE) 0, then the full (or any) FSM cannot be applied to any A € M{U, E, IE). We however 
show that shifting the system Ax = b down by k rows leads to a system to which the full (and 
hence any) FSM applies. 

In Section 4 we bound and compare the norms, and the norms of inverses, of bi-infinite, semi¬ 
infinite and finite Jacobi matrices over ([/, E, IE). This, in particular, allows the study of the 
resolvent norms, and hence the pseudospectra, of these operators and matrices. For example we 
show for A+ G 'I'E+([/, E, IE) and A G 'I'E(17, E, IE) that, analogously to the corresponding result 
for the spectrum, specPA+ = specf^ U G, for £ > 0 and p = 2. Here G is any of the sets discussed 
above that closes the gap between spec A = speCgg^ and spec . And we are able to make close 
connections between the pseudospectra of finite and infinite matrices, for example showing that the 
union of all finite matrix pseudospectra, Ujp^Mfi^(u,v,w)Spec§F, coincides with specPA+, for p = 2. 
Our results in this section are a substantial generalisation of results in a study [12] of spectra and 
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pseudospectra of a particular pseudoergodic Jacobi operator, the Feinberg-Zee random hopping 
matrix {U = W = {±1}, V = {0}). Our results on the relation between norms of inverses (and 
hence pseudospectra) of finite and infinite stochastic Toeplitz matrices, in particular our results on 
the convergence of norms, condition numbers, and pseudospectra of finite sections to their infinite 
stochastic Toeplitz counterparts, reproduce, in the case that U, V, and W are singletons, results 
for (classical) Toeplitz operators and matrices [51, 3, 5]. 

Main techniques. Besides the limit operator techniques behind Lemmas 1.1 and 1.2 that 
were the core of [39], our second main tool is a “glueing technique” - see (37) and (38) - that is 
used in the proofs of two of our main results. Theorems 2.2 and 3.3 as well as in their quantitative 
versions. Propositions 4.5 and 4.10. This latter technique often complements the earlier in terms 
of relating finite, semi- and bi-infinite matrices to each other. 


2 Spectra of Pseudoergodic Operators 

We recall, as discussed in Section 1, that throughout the paper we consider the operators in 
M{U,V,W) D 4'E(t/,F,W) and M+{U,V,W) D 4'E+(C/, F, W) as acting on £p(Z) and F(N), 
respectively, for some p G [l,c»], and that the invertibility, Fredholmness (and also the index if 
Fredholm) of these operators is independent of p. In Theorem 2.1 of [39] the following was shown: 


Proposition 2.1 The following statements are equivalent: 

(i) one operator in 'I'E([/, F, IF) is Fredholm, 

(ii) one operator in 'I'E([/, F, IF) is invertible, 

(Hi) all operators in M(U,V,W) are Fredholm, 

(iv) all operators in M(U,V,W) are invertible, 

(v) one operator in 'I'E_|_(17, F, IF) is Fredholm, 

(vi) all operators in M^(U,V,W) are Fredholm. 

All these equivalences follow quickly from Lemma 1.1. Since the occurrence of one (and hence 
all) of the properties {i)—{vi) is obviously not a matter of a concrete operator but rather of the 
interplay between U, V and IF, we will call the triple ([/, F, IF) compatible if {i)-{vi) hold. We 
will see below in Proposition 4.1 that, if (t/, F, IF) is compatible, then also the inverses in (iv) are 
uniformly bounded, in fact bounded above by llA“^|j for any A G 'I'E(17, F, IF). 

The equivalence of (i), (ii), {iv), (v) and (vi) can also be expressed as follows: for every 
A G 4'E([/, F, IF) and every B+ G 4'E+(t/, F, IF), it holds that 

speCe,,A = spec A = (J specC = [J speCe,sC+ = speCe,,S+. (13) 
ceM{u,v,w) c+eM+(u,v,w) 

Let us denote the set (13) by S([/, F, IF) since it clearly depends on U, V and IF only and not on 
the choice of A or B^. Then (U, V, IF) is compatible iff 0 ^ S(17, F, IF). 

To get a first idea of the set S := T,{U, V, IF), let us look at simple lower and upper bounds on 
E, our discussion here taken from [39, Theorem 2.1 a)]. A lower bound on 

E = (J specC (14) 

ceM(u,v,w) 

is clearly found by taking this union over a set of simple matrices C G M{U,V,W) for which 
spec C is known explicitly or is easily computed. Natural candidates are matrices C with periodic 
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diagonals, and the simplest among those are matrices with constant diagonals - so-called Laurent 
(or bi-infinite Toeplitz) matrices. 

If C is the (only) element of M({u}, {?;}, {w}) with some u G U, v G V and w G W, i.e. C is 
the tridiagonal Laurent matrix with u, v and w on its diagonals, then [4, 6] 

specC = {ut + v + wt~^ : t G T} =: E{u,v,w) (15) 

is the ellipse depicted in Figure 2.1. a below. Note that 


E{u,v,w) = V + E{u,w) with E{u,w) := E(u,0,w). 

Also note that (15) gives the ellipse E{u,v,w) an orientation, based on the counter-clockwise 
orientation of the unit circle T: the ellipse is oriented counter-clockwise if |u| > |r(;|, clockwise if 
|u| < jwl, and collapses into a line segment if |u| = |w|. From (14) and specC = E{u,v,w) if C is 
Laurent, we get that the union of all ellipses E{u, v, w) with u G U, v G V and w G W is a simple 
lower bound on S: 

y E{u,v,w) C 11(17,14, IF). (16) 

u£U,v£V^w£W 

Because we will come back to Laurent and Toeplitz operators, let us from now on write 

T{u,v,w) := uS + vl + wS~^ and T+{u,v,w) 

for the Laurent operator T G M({u}, {u}, {re}), acting on ^^(Z), and for its compression r+ to 
£P(N), which is a Toeplitz operator. Here we write S for the forward shift operator, S : x t-G y 
with y{j -I- 1) = x(j) for all j G Z, and 5'“^ for the backward shift. From (13) and (15) (or [4, 6]), 
speCgggT+(u,u,re) = specT{u,v,w) = E{u,v,w). Further, 

spec T+(u, u, w) = com E(u,v,w) (17) 


is the same ellipse but now filled [4, 6]. (Here convS' denotes the convex hull of a set S C C.) Let 
£’in(M, u, w) and i7out('w, v, w) denote the interior and exterior, respectively, of the ellipse E{u, v, w), 
with the understanding that Ein{u, v,w) =0 and i?out('w, v, w) = C\E{u, v, w) when |tt| = Iwl and 
the ellipse E{u,v,w) degenerates to a straight line. The reason why the spectrum of a Toeplitz 
operator T+ is obtained from the spectrum of the Laurent operator T (which is at the same time 
the essential spectrum of both r+ and T) by filling in the hole £’in(u,u,ry) can be found in the 
classical Coburn lemma [18]. We will carry that fact over to stochastic Toeplitz and Laurent 
operators. A key role will also be played by the following index formula. Let wind(r, z) denote the 
winding number (counter-clockwise) of an oriented closed curve F with respect to a point z ^ T. 
For 0 ^ E{u,v,w), so that T+ is Fredholm, it holds that [4, 6] 


indT+(M, u, w) = — wind(i7(u, u, ui), 0) 


0 , 0 € Eont{u,v,w), 

sign(|w| - juj), 0 G Ei^{u,v,w). 


(18) 


To get a simple upper bound on E, write A G 4'E([/, V, IF) as A = D + T with diagonal part 
D = diag(ui) and off-diagonal part T and think of A as a perturbation of H by T with ||T|| < e, 
where e •.= u* + w* and 

u* := maxlul, w* := maxlwl. 

«ec/ wcw 

Since A is in the e-neighbourhood of D, its spectrum spec A = S is in the e-neighbourhood of 
specZl C V. (Note that D is normal or look at Lemma 3.3 in [39].) In short, 

E(C/,F,VF) C V + {u*+w*)B (19) 
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(recall that D := {z £ C : < 1} is the open unit disk, and ID) is its closure). Note that the same 

argument, and hence the same upper bound, applies to the spectra of all (singly or bi-)infinite and 
all finite Jacobi matrices over 17, V and W. 

Sometimes equality holds in (19) but often it does not. For 17 = {1}, V = {0} and W = T, the 
lower (16) and upper bound (19) on S coincide so that equality holds in (19) saying that E = 2D. 
If we change W from T to { — 1,1} then the right-hand side of (19) remains at 2D while E is now 
smaller (it is properly contained in the square with corners ±2 and ±2i, see [12, 13]). Taking W 
even down to just {1}, the spectrum E clearly shrinks to [—2,2] with the right-hand side of (19) 
still at 2D. So the gap in (19) can be considerable, or nothing, or anything in between, really. 

Equality (13) contains the formula 

speCe,,S+ = IJ speCessC'+ = E 

C+&M+(U,V,W) 

for all 73+ G 'I'E+(17, V, W). One of our new results. Corollary 2.5 below, is that 

spec 73+ = spec (7+ =: E+ (20) 

C+&M+{U,V,W) 

holds independently of 73+ £ 'I'E+(7/, V, W). 

Upper and lower bounds on E+ = E+ (7/, V, W) can be derived in the same way as above for 
E. This time, because of (17), the ellipses in the lower bound (16) have to be filled in, while the 
upper bound from (19) remains the same, so that 

cony E{u,v,w) C Y,+ {U,V,W) C U-I-(u*-I- w*)D. (21) 

The results in this section will also make precise the difference between E and E+. 

For nonzero Toeplitz operators T+ (semi-infinite matrices with constant diagonals), acting 
boundedly on £^’(N), the following classical result fills the gap between essential spectrum and 
spectrum; at least one of the two integers, q:( 7+) and /3(T'+), is always zero. So if their difference 
is zero (i.e. T+ is Fredholm with index zero) then both numbers are zero (i.e. T+ is injective and 
surjective, hence invertible). This is Coburn’s Lemma [18], which was also found, some years ear¬ 
lier, by Gohberg [23] (but for the special case of Toeplitz operators with continuous symbol). Here 
is a new cousin of that more than 50 year old lemma: 

Theorem 2.2 If (U,V,W) is compatible (i.e. (i)-(vi) hold in Proposition 2.1) then every 73+ £ 
M+(t/, y,IU) is Fredholm and at least one of the non-negative integers q:(73+) and /3(73+) is zero. 


Proof. Let {U,V^W) be compatible and take 73+ £ M+(I/, y,IU) arbitrarily. Then 73+, with 
matrix representation (6^ is Fredholm since {i)-(vi) of Proposition 2.1 hold. Suppose that 

a(73+) > 0 and I3{B+) > 0. Then there exist x £ ^^(N) and y £ ^^(N), with p~^ -\- q~^ = 1, a: yf 0, 
and y ^ 0, such that 73+a: = 0 and 73jj/ = 0. Let a,b €C, define z £ i°°{Z) by 


z=(.. 


,ay2,ayi, 


0 , bxi,bx2,- ■ ■) 


T 


where the box marks the entry zq, and define 73 £ M{U, V,W)hy its matrix representation {bij)ij^z 
with 
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i,j £ N, 

l*-il > 1, 

i,j £ -N. 



The remaining entries 60 ,- 1 , 61,0 G U, 60,0 G V and 6 o,i, 6 _i,o € W of B may be arbitrary. Then 
B+x = 0 and B^y = 0 together imply that {Bz)j = 0 for all j ^ 0. Further, 

{Bz)o = bo-iayi + 60 , 16 x 1 . 

Clearly we can pick a, 6 G C with a 7 ^ 0 or 6 7 ^ 0 (so that z 7 ^ 0) to ensure that also {Bz)q = 0 
which implies that Bz = 0, so that B is not invertible on But this is a contradiction since 

B G M{U, V, W) is invertible by our assumption that the triple ([/, V, W) is compatible. ■ 

An immediate corollary of [39, Theorem 2.1] or Proposition 2.1 and our Theorem 2.2 is the fol¬ 
lowing Coburn lemma for pseudoergodic operators. This result reduces to the usual Coburn lemma 
for Toeplitz operators (at least to the special case in which the Toeplitz operator is tridiagonal) in 
the case that [/, V, and W are singleton sets. 

Corollary 2.3 //i?+ G 4'E+(17, C, IT) is Fredholm then at least one of a{B+) and I3{B+) is zero. 

Proof. If an operator in 'I'E+(C/, V, W) is Fredholm then {U, V, W) is compatible by Proposition 
2.1. Now apply Theorem 2.2 to B+. ■ 

Similarly to the situation for Toeplitz operators, one can now derive invertibility of operators 
in M+([/, y, W) from their Fredholmness and index. The additional result here that every B+ G 
M+([/, V, W) is Fredholm with the same index was first pointed out in [39], as a consequence of the 
main result from [46] (see our Lemma 1.1 e)), and will play an important role in our arguments. 

Corollary 2.4 Let {U,V,W) be compatible. Then every B+ G M+(17, E, IT) is Fredholm with the 
same index 

k(U,V,W) := indS+ G {-1,0,1}. 

Further, ifmAB+ = 0 then B+ is invertible; z/indi?+ = 1 then a{Bjf) = 1 and j3{B.y.) = 0; and 
i/indi?+ = —1 then a{Bj^) = 0 and (3{Bjf) = 1. These statements are independent of the choice 
of B+ G M+iU,V,W). 

Proof. If the triple {U,V,W) is compatible then every G M+{U,V,W) is Fredholm, by 
Proposition 2.1. To see that all B+ G M+{U, T, IT) have the same index, apply Lemma 1.1 e) to 
an A_|_ G '^Fj+{U,V,W). Because M+{U,V,W) includes the Toeplitz operators T+(u,v,w) with 
{u,v,w) G (17, y IT), it follows from (18) that k{U,V,W) G { — 1,0,1}. The remaining claims 
follow from Theorem 2.2. ■ 

Corollary 2.5 An operator 5+ G 'I'E+(t/, T, IT) is invertible iff all operators C+ G M+({7, T,1T) 
are invertible. In other words, for all B^ G \E'E_|_(17, V, IT), 

specil+ = specC+, 

c+eM+{u,v,w) 

this set denoted E+ = E+([/, T,1T) in (20). 

Proof. By Corollary 2.4 if one operator B+ G \E'E+(C/, T, IT) is invertible then (17,V,W) is 
compatible, k{U, V, IT) = 0, and all operators in M+{U, T, IT) are invertible. The formula for the 
spectrum follows by considering B+ — A/+ instead of B+. ■ 

As an extension to this corollary we will see in Theorem 3.1 that, if B+ G 'I'E+(t7, T, IT) is 
invertible, then 

sup llC}r^ll<oo. ( 22 ) 

c+eM+(u,v,w) 
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Before we explore further what the above results mean for spectra, we note some consequences 
of Corollary 2.4 that are essentially captured in [39, Theorem 2.4]. If (C/, V, W) is compatible then 
Corollary 2.4 tells us that every B+ G M+{U, V, W) has the same index k{U, V, W). In particular 
K = indr+(u,'u, w), for every {u,v,w) £ {U,V,W), this index given by (18). It follows that, if 
([/, V, W) is compatible, then this index must have the same value for all (u, u, w) £ (C/, V, W), so 
that either 0 G Eout{u, v, w) for all {u, v, w) £ {U, V, W), in which case /t = 0, or 0 G Ein{u, v, w) and 
jwj — juj has the same sign for all {u,v,w) G {U,V,W), in which case k = —wmd{E{u,v,w),0) = 
sign(|w| — |u|). Thus, if £ 4'E+(C/, V) W) is Fredholm but not invertible, either 

lu* := min IwI > u* (when k = 1) or u* := min juj > w* (when k = —1). (23) 

w^W u^U 

In the first case 0 is circumnavigated clockwise by all ellipses F(m, u, w), in the second the circum¬ 
navigation is counter-clockwise. 

For A G C, put F — A := {u — A : u G F} and note that B+ — A/+ G M+{U, F — A, IF) iff 
B+ £ M+{U,V,W). Similarly, B+-XI+ £ 4'F+(I/, F-A, IF) iff 5+ G 4'F+(I/, F, IF). We split the 
complex plane C into four pairwise disjoint parts. To this end, fix an arbitrary B+ £ 'I'F+([/, F, IF). 
The first part of the plane is our set E = S(f7, F, IF) from (13), 

E = {A G C : — A/+ is not Fredholm} 

= speCg,,B+ = y speCe,,C'+ 

c+eM+{u,v,w) 

= {A G C : (C/, F — A, IF) is not compatible}. 

The rest of the complex plane now splits into the following three parts: for k = —1,0,1, let 

Efe := {AgC : ind (B+- A/+) =/c} = {A G C : k(C/, F - A, W) = fc}. 

So we have a partition (i.e. a splitting into pairwise disjoint sets) of C: 

C = E U E_i U El U Eo, (24) 

s+ 


where the equality 

E+ = E U E_i U El = C \ Eo (25) 

holds by Corollary 2.4. So the difference between E+ = speci?+ from (20) and E = speCggg i?-|- 
from (13) is precisely E_i U Ei. 

The computation of these four parts of the plane, E, E_i, Ei and Eg, is of course far from 
trivial (otherwise spectral theory of random Jacobi operators would be easy) but we will compare 
this partition of C with another partition of C that is closely related and, in contrast, easy to 
compute. To do this, let £ denote the set of all ellipses E{u,v,w) with u £ U, v £ V and w £ IF, 
and put 
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Figure 2.1: a) The ellipse E{u,v,w) with general values u,v,w € C is derived from the zero-centered ellipse 
E{u,w) := E(u,0,w) after translation by v. The ellipse E(u,w) is centered at the origin and has orthogonal half¬ 
axes (the dotted lines) of length ||u.| di \'w\\, respectively, where the major axis bisects the angle between u and w at 
the origin. Depicted here is the ellipse E(u, w) for the particular values u = 3 and w = i. 

b) We see the splitting of C into the four parts E, -E_i, Eq and Ei for U = {—1,1}, V = {0} and W = {2}: the 
dark gray area is Ei, the light gray area is Eq and the rest (the white area plus the ellipse boundaries) is E. In this 
example, E—i is empty. Note that both ellipses, £'(—1,2) and £(1,2), are oriented clockwise. 


Eq ■= {XgC : 

X is outside of all ellipses in f}, 

En := {XgC : 

A is inside all ellipses in f}, 

El := (A G Efi 

: A is circumnavigated clockwise by all ellipses in £"}, 

E-i := (A G Efi 

: A is circumnavigated counter-clockwise by all ellipses in f}, 

E :— C \ {Eq U El U Fz—i), 

Eu := C\Eo = 

(J conv E{u,v,w). 


uGU,vGV,w^W 


\ 


(26) 


Obviously, at most one of the sets Ei and i?_i is nonempty. The set E consists of the points that 
lie on one of the ellipses, or they are inside some but outside other ellipses, or they are inside all 
ellipses but circumnavigated clockwise by some and counter-clockwise by others. We have 

C = E U E_i U ill U Eq, (27) 

in analogy to (24). But while the ingredients of (24) are in general notoriously difficult to compute, 
those of (27) are easily drawn. Before we relate (24) to (27), it is perhaps time for an example. 

Example 2.6 Take U = {—1,1}, V = {0} and W = {2}. Neither E nor E+ is precisely known 
in this case (but see [13, 27] for bounds on both). But the ingredients of (27) are easy to write 
down: draw all ellipses E(u,v,w) with u € U, v G V and w G W. In this case, there are only 
|17| • \V\ ■ \W\ = 2 ellipses: E(— 1,0,2) and E(1,0, 2). The situation is depicted in Figure 2.1.b 
(which is taken from [39]). The dark gray area is our set E^, the light gray area is Eq, and the 
rest (the white area plus the ellipses themselves) is E. Both ellipses are oriented clockwise since 
juj < jwj in both cases. So in this example, Ei = Ep, and i?_i is empty. □ 
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Now let us come to the relation between the partitions (24) and (27). From the discussion above 
(23) we see that 


S_i C i?-i, Si C i?i, So C Eq and hence i? C S. (28) 

So we have at least some simple upper bounds on S_i, Si and Sq and a lower bound on S. The 
upper bound on Sq is equivalent to the lower bound Ed on its complement S-i in (21). The lower 
bound FI on S is actually sharper than the lower bound (16). Further, from the discussion leading 
to (23) we see that 


f En, if w* > u*, 
[ 0 , otherwise, 


Ed, if > w*, 
0 , otherwise. 


(29) 


Recall that the difference between S_|. and S (i.e. the non-essential spectrum of 5+) is S_iUSi. 
From (28) we get that 


S^i :— S_i U Si C E—i U El —: F/±i C Ed- (^ 6 ) 

So the non-essential spectrum of B+ G 4>E_|_({7, V, W) is inside all ellipses in E. By (21), we have 

Ed G Ed G S+. (31) 

Now let us sum up. From (25), (30) and (31) we get 

S+ = SUS ±1 C SUEii C SUEn C SUEy C SUS+ = S+, 
so that all inclusions are in fact equalities and we have proven the following: 

Theorem 2.7 It holds that 

S+ = S U S ±1 = S U E±i = E U Ed = S U Ey (32) 

with Sj-i and E±i from (30) and Ed and Ed from (26). Moreover, exactly one of the following 

cases applies: 

i) tc* < u* and u* < w*, in which case S±i = E±i = 0 and Ed G S_|_ = S; 

ii) re* > u*, in which case S_i = E-i = 0 and Si C Ei = Ed; 

iii) u* > w*, in which case Si = Ei = 0 and S_i C i?_i = Ed- 

Of course Ed (and certainly Ed) is in general larger than the actual gap S+ \ S = S±i between 
the spectrum and essential spectrum, but equality (32) is still an attractive new bit of the picture: 
we do not know, very explicitly, what the sets S from (13) and S+ from (20) are, but we do now 
know explicitly what we have to add on to S to get S+. It has also recently been shown [27] that 
conv(i?u) = conY{E) is, surprisingly, both the closure of the numerical range and the convex hull 
of the spectrum for each operator in \E'E([/, V, W) U 'I'E_|_({7, V, W), and is hence a (very explicit) 
upper bound on both S and S_|_. Combining this result with Theorem 2.7 and (28) we have that 

E C S C S+ = S U Ey C conv(i?y) = conv(E). (33) 

Example 2.8 — Bidiagonal case. In [37] the bidiagonal case was studied, that means U = {0} 

or W = {0}. Let us say U = {0}. Then all our ellipses 

E{u,v,w) = E{0,v,w) = u-l-rcT, u G U, v gV, w gW 
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are circles with clockwise orientation. So we have that 

Eyj = [J ("c + w*]Qi) and En = ('c + w*D) 

v^V vGV 

with w* = maxwew I'u^l and w* = niinu,gvv |w|. In [37] it was shown that 

S = \ and hence, by (32), we have 1]+ = E^- 

So in this case, the partitions (24) and (27) coincide: 

E = E = Eij \ Er\, I]_i = E_i = 0 , Ei = E’l = i?n, Eq = Eq = <C \ E\j. 

For A S S+ \ S = Si = £^n, each B+ € 'I'E+(?7, F — A, IF) is Fredholm with index 1; precisely, 
a{B+) = 1 and I3{B+) = 0. □ 

Example 2.9 Take U = {!}, V = {0}, and IF = {0,2}. In this case there are two ellipses: 
£'(1,0) = T, with a counter-clockwise orientation, and 

£(1,2) = {tj.+ iti-.trjti €M. and {tr/3)"^+ tj = 1}, (34) 

with a clockwise orientation. Case i) in Theorem 2.7 applies, so that £u C S+ = S. Further, 
£(1,0) C conv£(l,2), so that £u = conv£u and, by (33), S-|- = E = £u = conv£(l,2). □ 


3 The Finite Section Method 

From Theorems 2.8 and 2.9 of [39], or see Lemma 1.2 and the comments immediately below that 
lemma (or [55]), we know that the full FSM (with /„ = 1 and = n) applies to a semi-infinite 
Jacobi matrix iff the operator itself and the set of associated semi-infinite matrices C- in (9) 
are invertible. To each matrix C_ corresponds (see Remark 1.3) a reflected matrix C+ = RCZR € 
M+{U, F, IF). Further, the set of all these reflected matrices C+ is all of M+{U, F, IF) iff A+ is 
pseudoergodic, as a simple consequence of Lemma l.ld). Similarly, the full FSM (with = —n and 
fn = n) applies to a bi-infinite Jacobi matrix A iff the operator itself and both sets of associated 
semi-infinite matrices £+ and C- in (9) are invertible and, again by Lemma l.ld), the union of 
the set of all matrices £+ and C+ = RCJR is the whole of M+{U, V, IF) iff A is pseudoergodic. 
As a simple consequence of these facts, and the results in sections 1 and 2, we obtain: 

Theorem 3.1 The following are equivalent: 

(a) the full FSM applies to one operator in d>E(U,V,W); 

(b) the full FSM applies to one operator in \1'E_|_([/, F, IF); 

(e) all operators in M+([/, F, IF) are invertible; 

(d) all operators in M^{U,V,W) are invertible, and the inverses are uniformly bounded; 

(e) one operator in 'I'E+(C/, F, IF) is invertible; 

(f) 0 ^ S+; 

(g) one operator in 'I'E+(C/, F, IF) is Fredholm with index 0; 

(h) all operators in M^(U,V,W) are Fredholm with index 0; 

(i) ([/, F, IF) is compatible and k(U, F, IF) = 0; 

(j) the full FSM applies to all operators in MlJJ, V, IF); 

(k) the full FSM applies to all operators in M+({7, F, IF). 
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Proof. The equivalence of (c), (g), (h) and (i) follows from Proposition 2.1 and Corollary 2.4. 
The equivalence of (a)-(c) is then clear from the remarks preceding this theorem and Lemma 1.2, 
as is the equivalence of (c) with (j) and (k) (or see Theorems 2.8 and 2.9 in [39]). That (c), (e) and 
(f) are equivalent is Corollary 2.5. To see the equivalence of (c) and (d) suppose that (c) holds, in 
which case also (j) holds. Then, by Lemma 1.2, all the operators 5+ and C_ in (9) arising from 
the full FSM applied to any A G M{U, V, W) are invertible and uniformly bounded. Thus, for all 
P e [l,oo], 

sup IjB^^llp < oo and sup IjCr^Hp < oo, 

so that also, where C+ = RCZR and since £'?(—N) is the dual space of .^p(—N) if p~^ + q~^ = 1, 
sup [jCliI^llq < oo. Since, see the remarks before the lemma, the collection of all operators 5+ and 
C+ is the whole of M+{U,V,W) in the case that A G 4'E(t/, V, IT), we see that we have shown 
(d), precisely that 

M+.p-.= sup ||B+^||p<oo, (35) 

B+€M{U,V,W) 

for 1 < p < oo. ■ 

Remark 3.2 a) So for pseudoergodic semi-infinite matrices A^ (the so-called “stochastic Toeplitz 
operators” from [57]), we get that the full FSM applies as soon as the operator is invertible. This is 
of course the best possible result since invertibility of A^ is a minimal requirement (it is necessary) 
for the applicability of the FSM. For (classical) banded Toeplitz operators T+, the same is true 
as was first shown in [24]. So in a sense, we also rediscover that classical result for tridiagonal 
Toeplitz matrices (by applying our result to the case when U, V and W are singletons). 

b) There is a similar coincidence between the FSM for pseudoergodic bi-infinite matrices (called 
“stochastic Laurent operators” in [57]) and for usual Laurent operators: in both cases, the FSM 
applies iff the operator is invertible and the corresponding semi-infinite principal submatrix (the 
Toeplitz part) has index zero. If the latter index is not zero, there is something that can be done. 
It is called “index cancellation”, and we will get to this in short course. 

c) Recall from Lemma 1.4 (also see [39, Thms. 2.8 & 2.9]) that the FSM applies with arbitrary 
monotonic cut-off sequences In and r„ if the full FSM is applicable. □ 

Recall that applicability of the FSM means that the truncated (finite) systems (7) are uniquely 
solvable for all sufficiently large n, say n > no, with their solutions Xn approximating the unique 
solution X of (6). The limit operator techniques behind Lemma 1.2 don’t reveal much about that - 
practically very relevant - number no; but in our pseudoergodic setting we can prove that actually 

no = 1 holds in (a), (5), (j) and (fc) of Theorem 3.1: 

Theorem 3.3 From the equivalent conditions (a)-(k) of Theorem 3.1 it follows that all finite 
Jacobi matrices over U, V and W, that means all F G Mf^n{U, V, W), are invertible. 

Before we come to the proof, let us note that this theorem implies 

spec J C E+ (36) 

for every Jacobi matrix J (finite or infinite) over U, V and VF, i.e., for every J G Mfin{U, V, W) U 
M+(C/, y, W) U M{U, V, W). Equality holds if J e TE+(C/, V, W). 

We now prepare the proof of Theorem 3.3. It combines a technique from the proof of [12, 
Theorem 4.1] with elements of the proof of Theorem 2.2. Let n G N. Given an n x n matrix 
F G M(in((J, V, W) and arbitrary elements u G U, v G V and w G W, we make the following 
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construction. Put 

/ •. 


B := 



w 


V 


e M{U,V,W), (37) 




where 0 marks the entry of B at position (0,0). We denote the semi-infinite blocks above and 
below 0 by B- and B+, respectively, so that 


/ 


B = 




\ 


w 


u w 

u 


B. 


(38) 


V y 

Precisely, with B = {bij)ij^z, we put B+ := (&y)ijgN £ M+{U, V, W) and B_ := (&y )ijg_N. 
Now, for a vector a; S C" and a complex sequence (rfe)fegz, put 


f : \ 


T-iX 


rox 


rix 


V 


/ 


leading to Bx = 


y 


\ 


- 2-1 


t-iFx 


Zo 


TqFx 


Zi 


tiFx 


Z2 


(39) 


where 0 and Zq mark the respective 0 positions and 

Zk = Tk^lUXn+TkWXi, 


k €Z. 


(40) 
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Lemma 3.4 For arbitrary Xi,Xn,u,w € C, there exists a sequence {rk)kei, in C so that, for the 
sequence [zk)k& from (40), either 

(r-o,r-i,...) e £°°(N) \ {0} and zi= Z 2 = ■■■ = 0 (41) 

or 

(...,r_ 2 ,r_i) e £°°(-N) \ {0} and z_i = z _2 = ... = 0. (42) 


Proof. The choice of (r^) differs, depending on whether (and which) parameters are zero. If none 
of xi,Xn, u, w is zero, put p := —uxn!(uixi) and := for fc G Z, so that Zfc = 0 for fc G Z and 
(41) holds if IpI < 1, (42) if \p\ > 1. If uXn = 0, put tq := 1 and := 0 for fc 7 ^ 0, in which case 
(41) holds. If wxi = 0, put r_i := 1 and rk := 0 for fc 7 ^ —1, in which case (42) holds. ■ 

Now we have all that we need: 

Proof of Theorem 3.3. Let F G Mn{U,V,W) for some n G N. Suppose F is singular. Then 
there is an a; G C" \ {0} with Fx = 0. Choose a sequence (r^) in C so that (41) or (42) holds. 
Then, in the notations of (38) and (39), either i?+x+ = 0 with x+ := i|N G f°“(N) or B-X- = 0 
with i_ := ij-N G ^°°(—N). In either case, there is a non-invertible operator (either or the 
reflection RBfR of i?_) in M+{U, V, W), which contradicts our assumption. So F is invertible. ■ 


So, in the case k{U,V, W) = 0, the full FSM (and hence the FSM with any monotonic cut-off 
sequences and r„) applies to every A G M{U,V,W) and to every A+ G M+{U,V,W), where 
unique solvability of the finite systems (7) already starts at n = 1. 

In the remaining cases (assuming, of course, that {U,V,W) is compatible, i.e. k{U,V,W) is 
defined), it holds that 

K := k{U,V,W) = ±1, 

and the FSM cannot apply - no matter how the cut-offs are placed (e.g. [38, Prop 5.2]). In the 
semi-infinite case the operator A^ is not even invertible since its index equals k 7 ^ 0. In the bi¬ 
infinite case, the way out is to move the system Ax = 5 up or down by one row (i.e. to renumber 
the infinitely many equations in ( 6 ) by increasing or decreasing the row number i by 1). So instead 
of Ax = 6 , the equivalent system S'^Ax = S'^b is solved, where S is the bi-infinite forward shift 
introduced earlier. Passing from A to A := S'^A preserves invertibility and corrects the index of 
the semi-infinite principal submatrices from k to 0. Indeed, 

indA+ = ind (S"'A)+ = ind (S"')+ -I- ind A+ = —k + k = 0. (43) 


This shifting process is called index cancellation; for Laurent operators A it goes back to [24, 33], 
for much more general operators, see e.g. [38, 39, 41, 55]. 

We claim that, after index cancellation, the full FSM applies to every A G M{U, V, W) also in 
the cases /t = ±1. This can be seen as follows: 


Case 1: k = +1. Then all ellipses E{u,v,w) with u € U, v € V and w € W are oriented 
clockwise, so that 


|w| > |u| > 0, u G U, w € W. 


Ax = b, i.e. 



■■■ ■■■ ^ 


( '■ ] 


( ] 

■ . V-2 W-2 


x{-2) 


b{-2) 

U-i V-i W-i 


x{-l) 


b{-l) 

0 

0 

0 


x{0) 

= 

6 ( 0 ) 

Ui Vi Wi 


xil) 


6 ( 1 ) 

U2 V2 


x{2) 


6 ( 2 ) 

■■■ j 


\ J 


1 ; / 


(44) 


(45) 
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to the equivalent system SAx = Sh, i.e. 


\ 

• ■ W-3 

V-2 W-2 

U-i V-i W-i 


f ] 

x(-2) 

x(-l) 

x(0) 


f ^ ^ 

b(-3) 

b(-2) 

b(-l) 

UO Vo Wo 


x(l) 


6 (0) 

Ui VI WI 


x(2) 


6 (1) 

■■■ ■■■ ■■■ } 


[ : J 


1 : / 


(46) 


V 

To see that the full FSM applies to the shifted system (46), it is sufficient, by Theorem 2.8 of [39], 
to show that all semi-infinite matrices of the form 


B~ = 


I 


\ 


Wi 



V2 

W2 


U3 

V3 

W3 


U4 

Vi 


W4, 


\ 


/ 


(47) 


with Ui G U, Vi G V and Wi G W are invertible on ^^’(N). So let 5+ be one of them. We start with 
injectivity: from ua ^ 0 for all z, by (44), we get successively a;(l) = 0, x{2) = 0,... as the only 
solution of Bj^x+ = 0. By [46], indi?+ = ind(S"‘A)+. But the latter is zero by (43), so that B+ 
is also surjective and hence invertible. 

Case 2: k = —1. Now all ellipses E{u,v,w) with u G U, v G V and w G W are oriented 
counter-clockwise, so that 

ju] > Jw] > 0, u G U, w G W. (48) 

Now pass from Ax = b, i.e. (45), to S~^Ax = S~^b, i.e. 

\ 


■■■ ■■■ 1 


1 1 


! : 

U-i V-i W-i 


x{-2) 


6 (-l) 

o 

o 

o 


x{-l) 


6 (0) 

Ui Vi Wi 


x{0) 

= 

6 (1) 

U2 V2 


x{l) 


6 (2) 

U3 


x{2) 


6(3) 



\ ) 


1 : 


(49) 


and check in a similar way, now using (48) and (43), that all semi-infinite matrices of the form 

\ 


B. = 


Ui Vi 

Wi 


U2 

V2 

W2 


U3 

V3 



Ui 


■■■ 


with Ui G U, Vi G V and Wi G W are invertible on £p(N). 

So also in case k{U,V,W) = ±1, the FSM applies, with arbitrary cut-off sequences In and r„, 
to every A G M{U, V, W) - but only after index cancellation. From (44) and (48) it is clear that 
every finite principal submatrix of the shifted (and therefore triangular) matrices in (46) and (49), 
respectively, is invertible. So, again, the finite systems (7) are uniquely solvable for all n (and not 
just for all sufficiently large n). 
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4 Norms, Norms of Inverses, and Pseudospectra 


In this section we bound and compare the norms and the norms of inverses of bi-infinite, semi¬ 
infinite and hnite Jacobi matrices over (U,V,W). The results are then expressed in terms of 
pseudospectra. 


4.1 Bi-infinite matrices 


We start with the simplest case: bi-infinite matrices in M{U, V, W). Not surprisingly, a prominent 
role is played by those in 4'E(t/, V, W). 

Proposition 4.1 Let A € 4'E([/, 1/,1V) and B G M{U,V,W) be arbitrary. Then ||yl|| > ||iJ||. If 
moreover {U,V,W) is compatible, then also ||4“^|| > 

Proof. First note that A and B are invertible by Proposition 2.1 if {U, V, W) is compatible. Now 
use Lemma 1.1 d), b) and a) - in this order. A £ \1'E([/, V, LE) implies B G a°P(A) by d), which 
then implies B~^ G a°'^{A~^) by b), so that ||i?|| < ||4|| and ||i?~^|| < ||4“^|| follow by a). ■ 


Corollary 4.2 For all A G 'I'E(t/, V, VE), we have 


Pll 


max ||i?|| =: M. 

BeM{U,V,W) 


If moreover (U, E, IE) is compatible then 


P-'II 


max \\B ^11 =■. J\f 
B&M(U,V,W) 


(50) 


If we have a particular p £ [1, oo] in mind, or want to emphasise the dependence on p, we will write 
Mp and Afp for the expressions M and Af defined in Corollary 4.2 (cf. (35)). 

The following proposition is a simple consequence of the observations that, ii A,B G BDO(fP(Z)) 
for all 1 < p < oo, and A = RB^R, where R is the reflection operator defined in Remark 1.3, 
then: (i) ||4||p = ||.B^||p = \\B\\q, for 1 < p < oo, if p~^ + = 1; (h) A G M{U,V,W) 

iS B G M{U,V,W)\ (iii) A is invertible iff B is invertible, and if they are both invertible then 

P-^llp=ll(S^)-'llp=l|i3-'ll,- 

Proposition 4.3 Forp,q G [l,oo], with p~^ + = 1, we have 

AAp = A4g and Afp = Afq. 

Proof. It is clear from the above observations that A4p = y 11^^ lip = ■Mq- Similarly 

Afp=Afg.m 


4.2 The relationship between semi- and bi-infinite matrices 


The semi-infinite case My{U, E, LE) is a bit more involved. We start with a simple observation: 

Proposition 4.4 a) For 4+ £ 'I'E+([/, E, VE) it holds that 

||4+|| = max ||R+|| = A4. 

B+eM+{U,V,W) 


b) If A+ G 'I'E+([/, E, VE) is invertible, i.e. {U,V,W) is compatible and n{U,V,W) = 0, then 

P;'II>AA. 
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Proof. Let A G M(U,V,W) be a bi-infinite extension of G P, fP). Then A G 

4'E(C/, P, fL), so that ||A|| = A4 and ||^“^|| = Af, by (50). Now ||A|| > ||A+|| since A+ is a 
compression of A, and ||A|| < ||A+|| since A G a°^{A+) = M{U,V,W). If A^ is invertible, this 
last fact implies that A~^ G so that ||^“^|| < ||A^^||. Finally, let B G M{U,V,W) be 

a bi-infinite extension of a given arbitrary G M+(?7, V, W). Then B G (t°p(A+) = M(U, V, W) 
and B^ is a compression of H, so that ||i3+ll — I|i3|| < P+ll- ■ 

A question we have been unable to resolve in general is whether ||Al|I^|| is the same for all A+ G 
tE'E+(C/, V, IE), and whether it is larger than M. We will see below that the following proposition, 
which is a partial complement of the bound ||A^^|| > M for A+ G 4'E+(17, E, IE), answers this 
question at least in the case p = 2. The arguments to obtain this proposition are a quantitative 
version of the proof of Theorem 2.2. 


Proposition 4.5 a) If A+,B+ G M+{U,V,W) are invertible and {U,V,W) is compatible, then 

min(||A+^||p, ||B+i||,) < Afp = Afq (51) 

holds for all p,q G [1, oo] with p~^ + q~^ = 1. 

b) Thus, if ([/, V, IE) is compatible with k(U, V, IE) = 0, so that all operators in V, W) 

are invertible, then (51) holds for all A_|_, B^ G M^{U, V, IE). 

c) //A+,i?_|_ G 'l'E+({7, E, IE) are invertible then equality holds in (51). 


Proof, a) If (C/, E, IE) is compatible then all operators in M{U, V, IE) are invertible and Afp = Afq 
holds by Proposition 4.3. Let A+,i?+ G M+(t/, E, IE) be invertible, and let C- := RB^R be the 
reflection of B+ as discussed in Remark 1.3. Abbreviate 

||A;1 ||p =: a and ||CI^||, = ||(RRXi?)-^||, = \\R{Bl)-^R\\p = ||(i?;i)^||p = ||i3;^||, =: b. 

Given an arbitrary e > 0, choose x G ^^’(N) and y G iP{—N) so that ||A+a;||p = 1, HC-j/Hp = 1 
with ||x||^ > aP — e and \\y\\^ > bP — e in case p < oo, and ||x||oo > a — e and ||?/||oo > 5 — e in case 
p = oo. Now put 



( 

C- 

w 

\ 


/ 

sy 

\ 

B := 

U 

H 

w 

G M ([/, E, IE) and x := 


0 




u 







\ 


A+ 


\ 

rx 

/ 


e F(Z), 


(52) 


where u G U, v G V and w G W are chosen arbitrarily, marks the entry at (0, 0) in B, J)j is 
at position 0 in x, and r,s G C are chosen such that {Bx)q = usy-i + wrxi equals zero, while r, s 
are not both zero. Then a? 7 ^ 0 and Bx = {sC-y, 0, rA+x)^ 7 ^ 0. Now, for p < 00 , 


, \\B-^Bx\\P imp \r\P\\j:\\P + \s\P\\y\\P 

||i? 2 ?||^ ||i? 2 ?||^ |r|P||A+:r||^ + |s|P||C'_ 2 /||^ 

= : taP + (1 — t)bP — e > min(a, b)P — e 


r\PaP + \s\PbP 
\r\P + |s|P 


holds, where t := |r|P/(|r|P + Isl^*) G [0,1]. Since this is the case for every e > 0, we conclude 
l|S-'llp > min(a, b). The case p = 00 is similar: 


l|S-^||oo > 
> 


||R~^Rx||oo _ Plloo _ max(|r|||a;||oc,, |s|||y||oo) 

||RS||oo ll^^?l|oo max(|r|||A+a;||oo, |s|||C'_y||oo) 

maxdrKa — e), |s|(6 — e)) max(|r|a, |s|6) . , 

-n I I i\- “ -/I I I IN-£ ^ min(o 

max(|r|, |s|) max(|r|, |s|) 
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Since e > 0 is arbitrary, again \\B ^||oo > min(a,6) follows. 

b) If, in addition, k{U, V, W) = 0 then all operators in M+{U, V, W) are invertible, so that a) 
can be applied to arbitrary A+, B+ G M+{U, V,W). 

c) If A+, G 4'E+([/, V, W) are invertible then the conditions of b) are satisfied, by Theorem 
3.1. In addition to b), note that, by Proposition 4.4 b), the minimum in (51) is greater than or 
equal to min(A/’p, A/",) if A+, B^ G 'I'E+(C/, V, W), which equals Mp = Nq, by Proposition 4.3. ■ 

We have seen already in Theorem 3.1 that A/+,p, our notation (35) for the supremum of HC^^Hp 
over all C+ G M{U, V, W), is finite for 1 < p < oo if (f/, V, W) is compatible with k{U, V, W) = 0. 
Propositions 4.4 b) and 4.5 imply that, in many cases, this supremum is in fact a maximum and 
coincides with the maximum (50) and with the norm of when A+ G 4>E+(I7, V, W). 

Proposition 4.6 Suppose that {U, V, W) is compatible with k{U, V, W) = 0, equivalently that all 
operators in M^(U, V, W) are invertible. Then, for every p,q G [1, oo] with p~^ + q~^ = 1; 

a) 

||A7^||p = max IICT^IL = A/+p = Afp for all A+ G (53) 
c+eM+{u,v,w) F 

or (53) holds withp replaced by q. If (53) holds we say thatp is favourable (for the triple {U, V, W)). 

b) p and q are both favourable iff M+^p = Af+,q. If P and q are both favourable, then 

||A;i||p = A4 ,p = Mp = Uq = M+,q = for all A+ G^E+{U,V,W). (54) 

In particular this holds for p = q = 2. 

Proof, a) Either ||7ll|l^|lp < Mp for all A^ G M+{U,V,W), or ||Al|l^||p > Mp for some A+ G 
M+{U,V,W). In the first case (53) follows immediately from Proposition 4.4 b). In the second 
case, by Proposition 4.5 b), < Mq for all i?+ G M+{U, V, W), and then (53), withp replaced 

by q, follows from Proposition 4.4 b). 

b) is an immediate corollary of a) and Proposition 4.4 b). ■ 

It is unclear to us whether every p G [l,oo] is favourable for every triple {U,V,W). Indeed, 
while, for every triple {U, V, W), p G [1, oo], and A+ G 4>E_|_([/, V, W), it follows from Propositions 
4.4 b) and 4.5 c) that 

ATp < ||A;i||p < M+,p, (55) 

it is unclear to us whether or not there are examples for which (55) holds with one or both “<” 
replaced by “<”. Likewise, it is unclear to us whether or not there are cases where ||Al|l^||p ^ 
||f?+^||p with A+,B+ G '^E+{U,V,W). A simple case to study is that where U, V, and W are 
singletons, in which case M+{U,V,W) = 4>E-|_(C/, P, W) has only one element A+, which is a 
tridiagonal Toeplitz operator. If A_|_ is invertible, then A/+,p = ||Al|l^||p. Example 6.6 of [4] shows 
a banded Toeplitz operator A+, for which ||Al|l^||p and ||Al|l^||g differ: if this were a tridiagonal 
banded Toeplitz operator then this would provide an example of a triple {U, V, W) where p and 
q are not both favourable. On the other hand, the following result shows that all p G [l,oo] are 
favourable for some classes of triples {U,V,W). 

Proposition 4.7 a) Suppose that p,q € [l,oo] and p~^ + q~^ = 1. Then p is favourable for 
{U, V, W) iff q is favourable for {W, V, U). 

b) If U = W or 0 G U U W, then, for all p G [1, oo], p is favourable (for {U, V, W)) and (54) 
holds. 

Proof, a) This is clear from Proposition 4.6, since A+ is in M+{U, V, W) or 'kE+([/, V, W) iff Aj_ 
is in M+(W, V, U) or 4'E+(W, V, U), respectively, and A+ is invertible iff Aj_ is invertible, in which 
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case P;^||,= ||(A|)-i|U. 

b) In the case that U = W \t follows from a) that p and q are both favourable. To show that they 
are both favourable if 0 G W, suppose that one of the two (say p) is not favourable. Then there 
exists A+ G M+{U, V, W) and x G £p(N) with ||yl+a;||p < A/^^||x||p. Take any B+ G M+{U, V, W), 
set C- = RBJR (as in the proof of Proposition 4.5), and define B G M{U,V,W) and x G ^^(Z) 
by (52), but choosing in particular w = 0 G W, r = 1, and s = 0. Then ||i?i||p = ||A+a:||p < 
A/'p'^lkllp = A/p-i ||i||p. But this implies that \\B ^||p > Afp, a contradiction. Thus p and q are 
both favourable if 0 G W, and it follows from part a) that they are both favourable also when 
0 G t/. ■ 

As one example, part b) of this proposition applies to the Feinberg-Zee random hopping matrix, 
A+ G 4'E+(C/, V, IF) or A G ^E(U, V, IF) with U = W = {±1} and V = {0} (e.g. [11, 13, 12, 27, 
28, 29]), so that (54) holds in that case (cf. [12, Theorem 3.6]). 

4.3 The relationship between finite and infinite matrices 

In this subsection we obtain more quantitative versions of the results of Section 3. We first note 
the following finite version of Proposition 4.3, proved in the same way, using the observation that, 
for every B G Mn{U,V,W), A = RnB^ Rn G M{U,V,W), where i?„ = (r^ is the nxn 

matrix with = Si^n+i-j, where Sij is the Kronecker delta. 

Lemma 4.8 For p,q G [l,oo] with p~^ + q~^ = 1 and n G N, we have M.n,p = A4„,g and 
Afji^p = SO that At fin,p — At fin,q and Affin^p — Afein^q! where 

Adn^p .— sup JjFjJp, Afn^p .— sup JJF' ][p, 

F&M„{U,V,W) FeM„{U,V,W) 

Atfin.p := sup l]Fllp and Affin,p ■= sup ||A“^l|p. 

FGMtiBU,V,W) F^Mnn{U,V,W) 


The following simple lemma relates Atfin,p to Atp, defined in Corollary 4.2. 

Lemma 4.9 For p G [1, oo], Atfin,p = lim„^oo At^^p = Atp. 

Proof. Let A G 4'E(?7, F, IF) so that UAjjp = A4p by Corollary 4.2. For F G M{in{U,V,W), 
llEjjp < |lA|]p, since every F is an arbitrarily small perturbation of a finite section of A. On 
the other hand, if An is the finite section of A given by (8), then we have noted in (10) that 
liminfn—^-oo ||An||p lidjjp. I 

The following is a more quantitative version of Theorem 3.3: 

Proposition 4.10 a) Properties (a)-(k) of Theorem 3.1 are equivalent to: 

(1) all F G Mfin{U,V,W) are invertible and their inverses are uniformly bounded. 

If (a)-(l) are satisfied then 

A/fin, p = max(A4,p , A/+,g), (56) 

for every p,q G [1, oo] with p~^ + q~^ = 1. 

b) In the case that p and q in a) are both favourable, (56) simplifies to 

A/fin,p — A/fin,(j = Af+,p — Af.\-^q = Afp = Afq. (57) 
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Proof, a) If (/) holds then, by the equivalence of ii) and hi) in Lemma 1.2 and the definition of 
stability, (e) holds. But this implies invertibility of all F € Mfin{U,V,W) by Theorem 3.3. The 
uniform boundedness of the inverses F~^ (and hence {1)) will follow if we can prove “<” in (56). 

To see that (56) holds, fixp G [1, oo], n G N, and an G M„([/, V, W). To estimate =: /, 

fix x G C" with ||fx||p = 1 and ||a;||p = /. As in the proof of Theorem 3.3, define B by (37) and 
i?_i_ and B_ as in (38), and define x by (39). Again choose (r^) as in Lemma 3.4 so that (41) or 
(42) holds. First assume it is (41). 

Case 1: p = oo. From (41) we get 

II \\B^^B+x+\\oo ||a+||^ supfcgKuW kfcllkllcx. / 1 

+ ||B+i+||oo ||B+x+||oo supfcgNulo} kfclll^ 2 ;||oo 1 ' 


Case 2: p < oo and x+ G £^’(N), i.e. {rk)t=Q £ ■^^(N). Then, by (41), 


\B- 


11 ^ > 


\\B+i+\\l 


ii^+iig sSkfenkiig fp 

WB+i+fp ES 


\\F-^ 


p 

p 


Case 3: p < oo and x+ ^ £^(N), i.e. (rfe)Eo ^ ^^(N). Then Sm '■= E^o —>■ oo as m ^ oo. 
Let m G N and put Xm '■= (^(l), 5(2), • • • , i((m + l)(n + 1)), 0,0, • • •) G fP(N). Then 




Er=oKi^ 


^F 


P 

m lip 


Efc=0 kfe|^l|5'2;||p + \urmXn\P Sm + \urmXn\P 


fP=\\F-^\\P 


since Sm ^ oo as m ^ oo and is bounded. 

So in either case we get HF’^^Hp < ||i3+^||p if (41) holds. The other case, (42), is analogous and 
leads to ||F'“^||p < ||i?Z^||p = HCfil^llg, where C+ := RB^R is the reflection of i?_ as discussed 
in Remark 1.3. Since we only know that (41) or (42) applies, but not which one of them, we 
conclude HF’^^Hp < max(||i?^^||p, HCfil^llg). Since F G M^n{U,V,W) is arbitrary and B+,C+ G 
M+{U,V,W), this finishes the proof of “<” in (56) and hence of the implication (a) => (/). The 
“>” in (56) follows from (10) or (11). 

b) follows from a), (53) and (54). ■ 

Combining Proposition 4.10 with Lemma 4.9 we can relate the p-condition numbers, condp(F') 
■= ll^’llp llT^^^llp of matrices F G Mfin{U,V,W), to the corresponding condition numbers of A G 
4'E(t/, y, IF) and A+ G 4'E+([/, V,W). For example, in the case that p and q are both favourable 
we have: 


Proposition 4.11 Suppose that {U,V,W) is compatible and k{U,V,W) = 0, so that every FSM 
given by (7) with In —>■ —oo and rn ^ oo is applicable to A € 4'E([/, F, IF), and every FSM given 
by (7) with Z„ = 1 and Vn ^ oo is applicable to A+ G \1'E_|_(C/, F, IF). Suppose also that p and q 
are both favourable. Then, for all these finite section methods, 

||A„||p —>■ ||A||p = ||A+||p = A4p, ||A„^||p —t ||A ^||p = ||A^^||p = J\fp, (58) 

so that condp(A„) —> condp(A) = condp(A+), 


as n ^ oo. Further, 

max condp(i3) = max condp(i3j_) = sup condp(F) = AlpAC, 

BgM{U,V,W) B+GM+{U,V,W) FGMnn{U,V,W) 

where the two maxima are attained, respectively, by all A G 'I'E(C/, F, IF), and by all A+ G 
4'E(t/, F, IF). 
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Proof. The first sentence is part of Theorem 3.1. And if p and q are both favourable, then the 
equalities in (58) follow from Corollary 4.2, Proposition 4.4 a), and (54). The stated limits are a 
consequence of 


Mp = ||A||p < liminf ||A„||p < limsup ||A„||p < Mfin,p = Mp 

n^oo n^oo 

and 

Np = Af+.p = 11 ^;^lip < liminf ||A“^||p < limsup ||A“^||p < = Np, 

n—>-oo n^oo 

by (10) and that, by Lemma 4.9 and Proposition 4.10 b), Alfin.p = Mp and A/fin,p = Mp. In 
the last displayed equation the first equality, and that these maxima are attained as stated and 
have the value MpAfp, follows from Corollary 4.2 and Propositions 4.4 a) and 4.6 b). That 
condp(F) < MpAfp is clear from A4fin,pA/fin,p = MpAfp', that in fact equality 
holds is clear from (58). ■ 

4.4 Pseudospectra 

We can rephrase our results on the norms of inverses, || J“^||, of Jacobi matrices J over {U, V, W) 
in terms of resolvent norms ||(J —A/)“^|| and pseudospectra, noting that J — XI is a Jacobi matrix 
over {U, V — W). In particular, J and J — XI are both pseudoergodic at the same time. In the 

language of pseudospectra. Corollary 4.2 and Proposition 4.3 can be rewritten as follows: 

Corollary 4.12 — bi-infinite matrices. For all A € 'I'E({7, V, W), e > 0 and p G [1, cx)], it holds 
that 

spec^A = := spec^i? and = E^, 

BeM{U,V,W) 

where p~^ + q~^ = 1. 

Summarizing Corollary 2.5, Theorem 2.7 and the results in Section 4.2, and recalling the nota¬ 
tions En and Eu from (26), we obtain: 

Proposition 4.13 semi- vs. bi-infinite matrices. 

a) For every A+ G \E'E_|_(17, V, W), e > 0 and all p G [1, oo], it holds that 

T.P UE+C specPA+ c ■= U specPC+. (59) 

c+eM+{u,v,w) 


b) For all A+, G M+([/, V, W), e > 0 and p,q € [1, oo] with p ^ + q ^ = 1, it holds that 

specPA+ n spec«.B+ c E^^^ n E« = Ef U G, (60) 

for each G G {if±i, An, Eu, E+}. Equality holds in (60) if G 'I'E+(t/, V, LE). If < u* 

and u* < w*, then E±i = 0 , so that (60) holds with G = 0 . 

c) //£ > 0 and p G [l,oo] is favourable, in particular if p = 2, then, for every A+ G 
^E+iU,V,W), 

specPA+ = EP_^ = EP U G, (61) 

for each G G {E±i, Ep, Ey, E+j. In particular (61) holds with G = 0 if w^: < u* and u* < w*. 

d) If p, q G [1, oo] with p~^ +q~^ = 1 are both favourable, in particular if Q & U AW or U = W, 
then (61) holds for both p and q, and spec§A_|_ = spec|A_|_. 
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Regarding the different possibilities for G in the above proposition, recall from Theorem 2.7 that 
each indicated choice closes the gap between E and E+ in the sense that E U G = E_|_, including 
the choice G = 0 if w* < it* and u* < 11 ;*. 

Proof. First, note that E C Ef since E = spec A C specfA = E^ for every A € 'I'E({7, V, W). 

a) Let A+ e «'E+(C/, V,W). If A e E+ = spec A+ then A £ spec§A+. If A G Ep \ E D E§ \ E+ 
then, for every A G 4'E(C7, E, IF), ||(A+ — A/+)“^||p > |1(A — A/)“^||p > s~^, by Proposition 4.4 

b), so that A G specPA+. 

b) We start with the equality in (60). By E C Ef and Theorem 2.7, EpuG = EpuEUG = 
EP U E+, for each G G {E±i, En, -Eu, E+}. 

To show that specPA+ n specfi?+ C Ef U E+, let A^,B+ G M+(U, V, IF) and A G specfA^ n 
spec|i?+, so that ||(A+ — A/+)“^||p > and ||(i3+ — A/+)“^||g > e~^. If one of these two 
numbers is infinite, i.e. A G specA+ or A G spec5+, then A G E+. If both ||(A+ — A/+)“^||p 
and 11(73+ — A/+)“^||g are finite and greater than £~^, consider A G 4'E(C7, F, IF), li A — \I 
is not invertible then A is in E C E+. li A — XI is invertible then, by (51), ||(A — A/)“^|| > 
min(||(yl+ — A/+)“^|lp, ||(i?+ - A/+)“^||g) > £“^, so that A is in Ef. 

Now suppose that A+,i3+ G \I'E+(17, F, IF) and that A G Ep U E+. If A G E+, then A G 
specA+ = spec 73+ by Corollary 2.5, so that A G specPA+ n spec®73+. If A G Ef then, by 
Proposition 4.5 c) and (50), min(||(Al+ — A7+)“^||p, ||(73+ — A7+)“^||q) = ||(A — A7)“^|| > 
for every A G 4iE([/, F, IF), so, again, A G specPA+ n spec®73+. In both cases it follows that 
AGE^^.nE^,,. 

c) This follows immediately from b) and (53). 

d) This follows from c), (54), and Proposition 4.7 b). ■ 

Here is the pseudospectral formulation of Proposition 4.10: 

Corollary 4.14 — finite, semi- and bi-infinite matrices. 

a) For all £ > 0 and p,q G [1, 00 ] with p~^ + q~^ = 1, it holds that 

SL+ := U U E®+,,. 

FGMti„{U,V,W) 


b) If p and q are both favourable, in particular if p = q = 2, this simplifies to 

= specPH+ = E^,, = EPUG = E®UG = E^,, = spec®73+, (62) 

which holds for all H+, 73+ G 'I'E+([7, F, IF), and for each G G {73±i, 73n, 73u, E+}, with E±i = 0 
if ic* < It* and it* <w*. 

Proof, a) This follows immediately from Proposition 4.10. In particular, 
sup ||(F-A7)-i||p 

FeM!in{U,V,W) 

= maxi sup ||(H+- A7+)"^||p , sup || (H+- A7+)“^ ||, 
\A+eM+{U,V,W) A+eM+{U,V,W) 

holds for all A G C, where both sides are infinite at the same time, 
b) follows from a) and Proposition 4.13 d). ■ 

Finally we study asymptotics of the pseudospectra of the finite sections of H+ G 4tE+(f/, F, IF), 
generalising results from [12] for the Feinberg-Zee random hopping matrix, and results of [5] which 
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apply when U, V, and W are singletons. For this purpose we briefly recall notions of set conver¬ 
gence. Let C®, C*" denote the sets of bounded and compact, respectively, non-empty subsets of C. 
For a € C and non-empty B C C, let dist(a, B) := infbgs \a — b\. For A,B G let 

dH{A,B) ■= max I sup dist(a, B), supdist(&, A) j . (63) 

VaeA bGB ) 

It is well known (e.g., [32, 26]) that •) is a metric on the so-called Hausdorff metric. For 
A,B£ C® it is clear that dniA, B) = dniA, B), so that dniA, B) = 0 iS A = B, and dff (•, •) is a 
pseudometric on C®. For a sequence {Sn) C and S S we write Sn^S if c?//(S'n,5') 0. 

This limit is in general not unique: if Sn S, then Sn^T iS S = T. 

A second related standard notion of set convergence is the following [32, Section 28], [26, Defini¬ 
tion 3.1]: for sets S C C, we write Sn ^ S ii lim inf Sn = lim sup = S', where lim inf Sn is the 
set of limits of sequences (Zn) C C such that Zn € Sn for each n, while lim sup Sn is the set of partial 
limits of such sequences. Both lim inf Sn and limsupS'„ are closed sets [26, Proposition 3.2], and 
it is clear (see [26, Proposition 3.5]) that, for every set sequence (S„), liminf Sn = liminf S„ and 
lim sup Sn = lim sup S'„. These two notions of convergence are very close. Precisely, for (5„) C C-®, 
liminf = lim sup Sn iff liminfS'„ = liminfS'„ = lim sup = lim sup and (see [32, p. 171] 
or [26, Proposition 3.6]), this holds iff Sn^S for some S S C*", in which case Sn —t S. Further, 
S' iff ^ S', and ^ ^ S' iff ^ S. Thus, for (S„) C C® and S € C^, S„ ^ S iff S„ ^ S. 

We introduce the following stronger notion of set convergence (cf., [32, p. 21]). Given S C C 
and a sequence of sets S„ C C, we will write S„ S if S„ C S for each n, and if every z G S is 
also in S„ for all sufficiently large n. In symbols Sn S means that S„ C S for each n and that 
S = UmeN n„>m Sn- The following lemma is immediate from this definition and the observations 
above: 

Lemma 4.15 If Sn C <C is a set sequence, S C C, and Sn S, then Sn S. If, additionally, 
each Sn is non-empty and S is bounded, then also Sn^S. 

In this last proposition we use again the notation (4). We remark that a weaker version of this 
proposition, namely that, under the same assumptions on p and q, limsup specf = SpeCgA_|_, 
can be deduced from general results on the finite section method for band-dominated operators 
[56], results which, as already noted, also lead to (11). 

Proposition 4.16 Let A+ G 4>E_|_(C/, V, W) and let (A„) denote its finite sections (8) with = 1 
and Tn —>■ -|-oo. If p,q G [1, oo] with p~^ -|- q~^ = 1 are both favourable, in particular if p = q = 2 
or OsC/UlF or U = W, then, for all e > 0, 

spec^An ^ spec£A+ = as n ^ oo, 

so that also spec^A^ specfA+ = SpeCgA+ and spec§A„-^spec^A_|_. 

Proof. By (62), specfA„ C spec§A_|_ = E]]. ^ for all n G N. Further, for every 0 < rj < e, 
it follows from [12, Theorem 4.4] that spec(]A_|_ C spec§A„ for all sufficiently large n, i.e. that, 
for some m G N, specPA+ c nn>mSpec§A„. Since Uo<j 7 <eSpec^A+ = specfA+, it follows that 
specf A+ = UmeN Iin>m spec§A„, and hence that specf A„ spec§A+ as n —5^ oo. The remaining 
results follow from Lemma 4.15. ■ 

Remark 4.17 the Hilbert space case p = 2 and the numerical range. We emphasise 
that our results are somewhat simpler in the important Hilbert space case p = 2. In particular, 
for A+ G 4'E+(17, V, W) and e > 0, where (A„) denotes the finite sections (8) of A+ with In = I 
and r„ -foo, we have from (62) and Proposition 4.16 that 

speCgA^ speCeA+ = = E^ U G, 
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for each G € {i?±i, E'n, £^u, S+}, with E±i = 0 if w* < m* and u* < w*. Further, in this Hilbert 
space case, as discussed above (33), Hagger [27] has shown that the closure of the numerical range of 
every operator in 'I'E(17, V, W) U 4'E-|_([/, V, W) is given explicitly by conv {E) = conv (Ey). Since, 
for a bounded operator on a Hilbert space, the £-neighbourhood of the numerical range contains 
the e-pseudospectrum (e.g. [57, Theorem 17.5]), and recalling (3), we can extend the lower and 
upper bounds (33) for the spectrum to give bounds on the 2-norm pseudospectra, that, for e > 0, 

Ey -b eD C C conv (Ey) -|- eD and E -b eD C C conv {E) -b eU. 

Example 2.9, with U = {1}, V = {0}, W = {0, 2} and Ey = convEy, is a case where these bounds 
on g are sharp: in that case = E^ ^ = conv (i?y) -beD = convi?(l, 2) -beD, with E(l, 2) the 
ellipse’ (34). □ 
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